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Abstract. We investigate the numerical approximation of (discontinuous) entropy 
solutions to nonlinear hyperbolic conservation laws posed on a Lorentzian manifold. 
Our main result establishes the convergence of monotone and first-order finite volume 
schemes for a large class of (space and time) triangulations. The proof relies on a dis- 
crete version of entropy inequalities and an entropy dissipation bound, which take into 
account the manifold geometry accurately and generalize techniques and estimates that 
were known in the (flat) Euclidian setting, only. The strong convergence of the scheme 
then is then a consequence of the well-posed theory recently developed by Ben-Artzi 
and LeFloch for conservation laws on manifolds. 



1. Introduction 

In the present paper, we consider discontinuous solutions to nonlinear hyperbolic con- 
servation laws posed on a globally hyperbolic Lorentzian manifold. Our main objective 
is to introduce a class of first-order and monotone finite volume schemes and derive 
geometrically natural and nonlinear stability properties satisfied by these schemes. In 
turn, we will conclude that the proposed finite volume schemes converge (in a strong 
topology) toward the entropy solutions recently defined in Ben-Artzi and LeFloch [2J, 
who established the well-posedness theory for conservation laws posed on (Riemannian 
or Lorentzian) manifolds. Our proof can be regarded as a generalization to Lorentzian 
manifolds of the method introduced in Cockburn, Coquel and LeFloch [3] for the (flat) 
Euclidean setting and extended to Riemannian manifolds in Amorim, Ben-Artzi, and 
LeFloch pQ. 

Major difficulties arise in working with a partial differential equation on a Lorentzian 
manifold, rather than on the customary (flat) Euclidian space. A space and time triangu- 
lations must be introduced and the geometry of the manifold must be taken into account 
accurately in the discretization. In addition, one cannot canonically choose a preferred 
time foliation in general, so that it is important for the discretization to be robust enough 
to allow for a large class of foliation and of space and time triangulations, that assume 
only limited regularity. From the numerical analysis standpoint, it is challenging to design 
and analyze schemes that are indeed consistent with the geometry of the given Lorentzian 
manifold. 

Our assumption of global hyperbolicity of the Lorentzian background allows us to en- 
compass physically realistic situation (arising, for instance, in the applications to general 
relativity). The class of schemes proposed in the present paper is quite general, and es- 
sentially assumes a monotonicity property on the flux functions only. Moreover, we over 
also a large class of space and time triangulations in which element can degenerate in 
certain directions of space or time. 
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Several new difficulties appear when trying to generalize to Lorentzian manifolds the 
convergence results in [TJ [3] . Our guide in deriving the necessary estimates was to ensure 
that all of our arguments are intrinsic in nature and (as the physics or geometry imposes 
it) do not rely explicitly on a choice of local coordinates. 

We show here that the proposed finite volume schemes can be expressed as a convex 
decomposition, and we derive a discrete version of entropy inequalities as well as sharp 
estimates on the entropy dissipation. Strong convergence towards an entropy solution 
follows from DiPerna's uniqueness theorem [3J. 

For further work on conservation laws on manifolds, we refer the reader to a pioneering 
paper by Panov [8], as well as to LeFloch and Okutmustur [7] and LeFloch, Neves, and 
Okutmustur [6]. 

An outline of this paper follows. In Section [2] we briefly state some results about 
conservation laws posed on manifolds. In Section [3] we present a formal derivation of the 
finite volume schemes and state our assumptions and main convergence result. Section 
H] is devoted to establishing stability estimates which are of independent interest but are 
also the basis to the convergence proof presented in Section [5j Finally, in Section [6l we 
present some examples of discretizations and particular schemes. 

2. Preliminaries 

2.1. Conservation law on a Lorentzian manifold. We will need some results on the 
well-posedness of hyperbolic conservation laws posed on Lorentzian manifolds., for which 
we refer to [2|. 

Let (M. d+1 ,g) be a time-oriented, (d + l)-dimensional Lorentzian manifold. Here, g is 
a metric with signature (—,+,...,+), and we denote by T p M rf+1 the tangent space at a 
point p G M d+1 . Recall that tangent vectors X (and more generally vector fields) on a 
Lorentzian manifold can be separated into time-like vectors (g(X,X) < 0), null vectors 
(g(X, X) = 0), and space-like vectors (g(X, X) > 0). The manifold is time-oriented, i.e. a 
consistent orientation can be chosen throughout the manifold, so that we can distinguish 
between past-oriented and future-oriented vectors. To the metric g one associates its 
Levi-Cevita connection (covariant derivative) V, which allows us to define the divergence 
div ff operator, as in the Riemannian the trace of the covariant derivative of a 

vector field. By duality arguments, the divergence of vector fields is known to depend 
only on the volume form associated with g. However, to assume a Lorentzian metric is 
most convenient to apply the terminology arising in general relativity when dealing with 
the initial- value problem (see below). 

Recall that, by definition, a flux on the manifold M rf+1 is a vector field f(u,p) G 
TpM 6 ^ 1 depending on a real parameter u. The conservation law on (M rf+1 , g) associated 
with a flux / is 

div 9 (f(u,p)) = 0, u:M d+1 ^K. (2.1) 
It is said to be geometry compatible if / satisfies the condition 

div ff /(«,p) = 0, uel,peM w . (2.2) 

Furthermore, / is said to be a time-like flux if 

g{d u f(u,p),d u f(u,p)) < 0, p G M d+ \ ueR. (2.3) 

We are interested in the initial- value problem associated with (|2.ip . We fix a space- 
like hypersurface !Ko C M d+1 and a measurable and bounded function uq defined on "Kq. 
Then, we search for a function u = u{p) G L°° (M. d+1 ) satisfying (|2.ip in the distributional 
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sense and such that the (weak) trace of u on "K coincides with it : 

n|^ = u . (2.4) 

It is natural to require that the vectors d u f(u,p), which determine the propagation of 
waves in solutions of (|2.ip . are time- like and future-oriented. 

We assume that the manifold M d+1 is globally hyperbolic, in the sense that there exists 
a foliation of M d+1 by space-like, compact, oriented hypersurfaces "Kt (i G R): 

M d+1 = \J ■Kt. 

tew. 

Any hypersurface "Kt is referred to as a Cauchy surface in M^" 1 " 1 , while the family "Kt 
(t G R) is called an admissible foliation associated with Jft . The future of the given 
hypersurface will be denoted by 

M d + 1 := |J -Kt. 

Finally we denote by n* the future-oriented, normal vector field to each "Kt, and by g l the 
induced metric. Finally, along "Kt, we denote by X 1 the normal component of a vector 
field X, thus X 1 := g^X,^). 

Definition 2.1. A flux F = F(u,p) is called a convex entropy flux associated with the 
conservation law (|2.ip if there exists a convex function U : R — > R such that 



F(u,p) = / d u U(u) d u f(u',p) du, p <E M d+1 , u£l. 

j4 measurable and bounded function u = u(p) is called an entropy solution of conservation 
law (|2.ip - (12. 2p if the following entropy inequality 



g{F{u),V g <f>)dV g + \ (div g F){u)<f>dV g 
M d+i Jm^ 1 



+ / go{F(u ),n ) <j>^ dV go - / U'(u){&w g f)(u)<j)dV g > 



3<o JM 1 ^ 1 

holds for all convex entropy flux F = F p (u) and all smooth functions (p > compactly 
supported in M^ +1 . 

In particular, the requirements in the above definition imply the inequality 

div 9 (F(u)) - (div g F)(u) + U'(u)(dw g /)(«) < 

in the distributional sense. 

2.2. Well-posed theory. We will use the following result (see [2]): 

Theorem 2.2. Consider a conservation law (|2.ip posed on a globally hyperbolic Lorentzian 
manifold M rf+1 with compact slices. Let "Kq be a Cauchy surface in M rf+1 , and uq : "Kq — ► 
R be a function in L°°(J{o). Then, the initial-value problem (|2.ip - (|2.4p admits a unique 
entropy solution u = u(p) G L^ c (Mi +1 ). Moreover, for every admissible foliation 'Kt, 
the trace u^ t G L 1 (!Kf) exists as a Lipschitz continuous function oft. When the flux if 
geometry compatible, the functions 

are non-increasing in time, for any convex entropy flux F. Moreover, given any two 
entropy solutions u, v, the function 

H/*0W - /*(w|M t )l|ii(Mt) 
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is non-increasing in time. 

Throughout the rest of this paper, we therefore assume that a globally hyperbolic 
Lorentzian manifold is given, and we tackle the problem of the discretization of the initial 
value problem associated with the conservation law (|2. 1 j) and a given initial condition 
where uq G L°°(CKo). We do not assume that the flux / is geometry compatible, and on 
the map div g / we require the following growth condition: there exist constants Ci, C2 > 
such that for all (u, p) G K x M 

|(div fl /)(«,p)| < C x + C 2 \u\. (2.5) 

Two remarks are in order. First of all, the terminology here differs from the one in 
the Riemannian (and Euclidean) cases, where the conservative variable is singled out. 
The class of conservation laws on a Riemannian manifold is recovered by taking M a!+1 = 
R x M, where M is a Riemannian manifold and f(u,p) = (u, /) G K X T p M, which leads 
to 

div g (/(it)) = d t u + div 9 (/(«)) . 

Second, in the Lorentzian case, no time-translation property is available in general, con- 
trary to the Riemannian case. Hence, no time-regularity is implied by the L 1 contraction 
property. 

3. Notation and main results 

3.1. Finite volume schemes. Before we can state our main result we need to intro- 
duce several notations and motivate the formulation of the finite volume schemes under 
consideration. 

We consider a (in both space and time) triangulation 7 h = IJifeT' 1 on com- 
posed of bounded space-time elements K satisfying the following assumptions: 

• The boundary dK of an element K is a piecewise smooth d-dimensional manifold 
without boundary, dK = Uecdft" e > wnere each cf-dimensional face e is a smooth 
manifold and is either everywhere time-like, null, or space-like. The outward unit 
normal to e G dK is denoted by n^^. 

• Each element K contains exactly two space-like faces, with disjoint interiors, 
denoted and e^, such that the outward unit normals to K, n R e + and n R e - , 
are future- and past-oriented, respectively. They will be called the inflow and the 
outflow faces, respectively. 

• For each element K, the set of the lateral faces d°K := dK\{e~^, e^} is nonempty 
and time-like. 

• For every pair of distinct elements K, K' G 7 h , the set KnK' is either a common 
face of K, K' or else a submanifold with dimension at most (d— 1). 

• "Ko C UxGT h ®K-> wnere ^-0 is the initial Cauchy hypersurface. 

• For every K G 7h, diame^- < h. 

Given an element K, we denote by K + (resp. K~) the unique element distinct from K 
sharing the face e\ (resp. e^) with K, and for each e° G d°K, we denote K e o the unique 
element sharing the face e° with K. 

The most natural way of introducing the finite volume method is to view the discrete 
solution as defined on the space-like faces separating two elements. So, to a particular 
element K we may associate two values, and associated to the unique outflow and 
inflow faces e]t, e~^. Then, one may determine that the value uk of the discrete solution 
on the element K, is the solution determined on the inflow face e~^ (one could just as 
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well say that uk is the solution uj- determined on the outflow face e~£, or some average 
of the two, as long as one does this coherently throughout the manifold). 

Thus, for any element K, integrate the equation (|2.ip . apply the divergence theorem 
and decompose the boundary dK into its parts e^-, e^, and d°K: 



9 P (f(u,p),n Ke +(p))dp- J _ g p (f(u,p),n Ke -(p))dp 

'\ ^ (3-1) 

+ iJ / 9 P (f(u,p),n K:e o(p))dp = 0. 

Note the minus signs on the first two terms. Indeed, for a Lorentzian manifold, the 
divergence theorem reads 



/ div g fdVn= / g{f,h)dV dn , 
Jn Jan 



with n is the outward normal if it is space-like, and the inward normal if it is time-like. 

Given an element K, we want to compute an approximation of the average of u(p) 
in the outflow face et-, being given the values of u^- on e^- and of for each e° E d°K. 

The following notation will be useful. Let / be a flux on the manifold M d+1 , K an 
element of the triangulation, and e C dK. Define the function {J K : R — > R by 

^K,e( U ) := /j^9p{f( u iP)i n K,e{P)) d P- ( 3 - 2 ) 

Also, if w : M d+1 — > R is a real function, we write 

Me :=jw{p)dp. 

Using this notation, the second term in (|3.ip is approximated by 

/ g P {f(u,p),n Ke -(p))dp~ \e K \^ f KeK (u K ), 
and the last term is approximated using 

g P (f{u,p),n e o(p))dp~\e°\ q K<e o {u~ K ,u K a ), 



where to each element K, and each face e° € d°K we associate a locally Lipschitz numer- 
ical flux function q^ ie o(n, v) : R 2 — > R satisfying the assumptions given in (|3.8p - ()3.10p 
below. 

Therefore, in view of the above approximation formulas we may write, as a discrete 
approximation of (|3,ip . 

\ e K\^ K+te +{ U K) ■= \ e K\V f K ,eS U K) ~ S l^l^.eofeXo)' ( 3 ' 3 ) 

which is the finite volume method of interest and, equivalently 

n i := (4+e+)" 1 (j^4, e -^)- E ^V(%.% ))' ( 3 - 4 ) 
v |e^| k e0< _ d0K \e K \ 

The second formula which may be carried out numerically (using for instance a Newton 
algorithm) is justified by the following observation: 

Lemma 3.1. For any K e 7 h , the function u i— > _ (u) is monotone increasing. 

K,e K 
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Proof. From f|3.2[) we deduce that 

duU*. (it) = I g P (d u f{u,p),n Ke - {p))dp > 0, 
since d u f(u,p) is future-oriented and n fc e - is past-oriented. □ 
If C !Ko , the initial condition (|2.4p gives 

o e K 

Finally, we define the function u h : M d+1 — > E by 

u h (p) := p E K. (3.6) 
For e G <9X we introduce the notation 

/ e (n,p) := g p {f{u,p),n K , e (p)). (3.7) 

3.2. Assumptions on the numerical flux- functions. We are now in a position to 
state our main result. First, we need to present some assumption on the flux function 
and the triangulations. 

The numerical flux functions q K e o(u,v) : M. 2 i— > R in the equation f|3.3|) satisfy the 
following assumptions: 

• Consistency property : 

<lK,eo(u,u) = ^ f e a{u,p)dp = H S Ke0 {u). (3.8) 

• Conservation property : 

<lK, e o(u,v) = -<lK s0 ,e° (v,V,), U, V G R. (3.9) 

• Monotonicity property : 

d u <lK,e a (u,v) > 0, d v q Kt£ o(u,v) < 0. (3.10) 

Remark 3.2. In fact, as we shall see in the examples in Section [6j it is natural, and even 
unavoidable, that the numerical flux-functions q^ e o depend on the whole geometry of 
the element K as well as of its neighbour, K e o . However, to keep the notations as simple 
as possible, we will use the above notations throughout. 

For each element K, define the local time increment 

\K\ 



\ e K\ 

We suppose that r := maxT/<- — > when h — > 0, and that for every K, 

h 2 /r -> 0. (3.11) 

For stability, we impose the following CFL condition which should hold for all K G 
7 h , e° G d K: 



\d°K\ 



sup 



d u ii f Ke0 (u) \ < (sup du{^ K+e+ r\u)) x <co, (3.12) 



K+ 
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3.3. Assumptions on the triangulations. We now describe the time-like geometric 
structure of the triangulation 7 . We will introduce a (global and geometric in nature) 
admissibility condition (see (|3.13p below) involving only the time evolution of the trian- 
gulation, which will be completely independent of its structure on space-like related faces. 
We stress that our method poses almost no restriction on the space-like structure of the 
discretization. 

We will define the local Cartesian deviation associated with a pair of elements K, 
K~, which is a quantity depending only on the geometry of each such pair of space- 
time elements, and propose an admissibility condition involving a global bound on this 
quantity. The local Cartesian deviation measures the amount by which the time evolution 
of the triangulation deviates from a uniform, Cartesian evolution, and is defined as follows. 
Let K G 7 h , and suppose first that its space-like outflow face e K satisfies the property 
that its center of mass, denoted by p K , lies on e K . Consider also the center of mass of 
d°K, denoted by p K . Then, define the bilinear form £(K) on (T p + M) 2 by 

where the vector w G T + M is the future-oriented unit tangent vector at p~t to the unique 

Pk " 

(if h is small enough) geodesic connecting p K to p K . So, if X, Y are two vectors defined 
at the point p K , we have, by definition, 

£ (K) (X, Y) = — w K <g> n K e+ (X, Y) = —g+ (X, w K ) g + (Y, n Re+ ), 

' k tk k k ' k 

and we define the local Cartesian deviation associated with K, K~ , as the form 

\K\E(K) - \K-\E(K~). 
This tensor-like operator measures the rate of change of the quantity 

\K\E(K)(X,Y) 

with respect to the time-like direction defined locally by the faces e K . Our admissibility 
criterion then states that this rate of change must not blow up faster than h~ l after 
summation over K G 7 h , as well as some flatness conditions on the space- like faces e K . 

Definition 3.3. We say that 7 h is an admissible triangulation if the following con- 
ditions are satisfied: 

(1) For every smooth vector field $ with compact support and every family of smooth 
vector fields ^k, K G 7 h , the Cartesian deviation satisfies 

V {\K\E{K)-\K-\E{K-)){$^ K ) < ^l||$|| L oosup||^|| i0 o, (3.13) 

— n K 

K"£T h 

for some function rj(h) such that 

r](h) -» 0. 

(2) For every K G 7 h , 

d g (p K ,e+)<\e+\T, (3.14) 

where d g denotes the distance function associated to the metric g. 

(3) For all smooth vector fields X defined on e K , 

\\g p (X{p),n K Mp))\\ c2{e + ) < \\X(p)\\ c2{e +y (3.15) 
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In the case where p^ et, the condition (|3.14|) means that as h — > 0, the center 
of mass p^ of approaches ej~-, even when we blow it up by a factor of order |e^|, 
as one would expect to find if e^- becomes flat in the limit. Under this assumption, 
we may uniquely extend the normal vector field + by geodesic transport to a small 

neighbourhood V of e^, in such a way that p~^ € V. Thus, in formula (|3.13j) . we can 
meaningfully speak of the vector n Ke +(p^) € T p +M. The condition (|3.15p is intended 

to rule out oscillations on the normal vector field due to the geometry of et. 

The assumption (13. 13ft amounts to a global geometric condition on the triangulation's 
Cartesian deviation. As we shall see in Section [6l this allows us to cover a large class of 
interesting, implementable space-time triangulations. 

3.4. Main result of convergence. Finally, we are in a position to state: 

Theorem 3.4 (Convergence of the finite volume method). Let u h be the sequence 
of functions generated by the finite volume method (|3.3h - (13.6P on an admissible triangu- 
lation, with initial data uq € L°°(3fo), and with numerical flux-functions satisfying the 
conditions (|3.8p - (|3.10p . and the CFL condition (|3.12p . Then, for every slice "Kt> ^ e se- 
quence u h is bounded in L°°( Use[o t] -^s), and converges almost everywhere when h — * 
towards the unique entropy solution u € L°°(M+ +1 ) of the Cauchy problem (12.11) . (|2.4p . 

In Section H] below, we will derive the key estimates required for a proof of Theorem 13.41 
which will be the content of Section [5j We follow here the strategy originally proposed by 
Cockburn, Coquel and LeFloch |3j. New estimates are required here to take into account 
the geometric effects and, especially, during the time evolution in the scheme. We will 
start with local (both in time and in space) entropy estimates, and next deduce a global- 
in-space entropy inequality. We will also establish the L°° stability of the scheme and, 
finally, the global (space-time) entropy inequality required for the convergence proof. 

Remark 3.5. Our formulas do coincide with the formulas already known in the Rie- 
mannian and Euclidean cases. In these cases, the function _ (u) is always the identity 

K,e K 

function. Therefore, the finite volume scheme reduces to the usual formulation found in 
[Ud]. Also, our expression for the time increment r and the CFL condition (|3.12p give 
the usual formulas when particularized to the Euclidean or Riemannian setting. 

Remark 3.6. More generally, and in the interest of practical implementation, one may 
replace the right-hand side of the equations (13.21) and (|3.5f) with more realistic averages. 
For instance, one could take an average of g(f(u,p), YiK,e) over iV spatial points pj given 
from some partition e 3 of e, 

1 N 

Ae( U ) = T-\^W\9{f{ u ^Vj)^KAPj))- 
161 3=1 

In view of these remarks, we see that the finite volume method may be given more 
generally by the algorithm which consists of fixing an averaging operator \i _ , and 

using the equation (|3.4p to iterate the method, with initial data given by 

, «(> 



u K := [i 



K 



However, any such average is just an approximation of the integral expression used in 
(|3.2p . This approximation can be chosen to be of very high order on the parameter h, 
by choosing appropriate quadrature formulas. Therefore, for the sake of clarity, we will 
present the proofs with \J K defined by (|3.2p and omit the (straightforward) treatment 
of the error term issuing from this approximation. 
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4. Discrete entropy estimates 

4.1. Local entropy dissipation and entropy inequalities. First of all, let us intro- 
duce some notations which will simplify the statement of the results as well as the proofs. 
Define 

With this notation, the finite volume method (|3.4p reads 

^kI^K^k) = \Zk\^k( u k) ~ Yl \ e °\ ( lK,e ( u K, u K e0 )- ( 4 -l) 

e°ed°K 

As in [3] and pQ, we rely on a convex decomposition of ^(lit-), which allows us to control 
the entropy dissipation. 

Define jx^ g0 by the identity 

\d°K\ 

#jr,e° := Vk(u k ) - — — (qK ie o(%,% ) ~<IkM u K> u k))i 



and define 



\et\ Jk 



Then, one has the following convex decomposition of //t(itt), whose proof is immediate 
from (|4.ip . 

^ {u+k) = Wk\ E l e °l^ e o- (4-3) 

Lemma 4.1. Let (U(u),F(u,p)) be a convex entropy pair (cf. Definition \2.1\) . For each 
K and each e = e^, e^, let VK, e '■ K — > R be the convex function defined by 

VkM ■■= ^A^ke)' 1 ^)), M G M- (4.4) 

Then there exists a family of numerical flux-functions Qk e (. u i v ) : — M satisfying the 
following conditions: 

• Qi<,e° i- s consistent with the entropy flux F: 

Q A>0 (u,u) = ^ e0 {u), K <E 7 h , e° € d°K,ue M. 

• Conservation property: 

QK, e o{u,v) = -Q Ke0>e o(v,u), u,v £R. 

• Discrete entropy inequality: 

\d°K\ _ _ (4.5) 

+ ^n-(Qv(%,% ) - QK, e o(u K ,u K )) < 0. 

\ e K\ 

At this juncture, we conclude from the inequality (14.51) that 

V K+,e+(T i K^ ~ V K+,e+(^K( U K)) 

\d°K\ _ _ (4-6) 

+ 1+1 {QK,e°( u Ki U K e0 ) ~ Qk,A U K^ u k)) ^ R K,e°> 
\ e K\ 
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where g0 is given by 

Proof. First of all, note that using (|4.4p we may write the inequality (|4.5p equivalently as 
F If. J \-l/r,+ ^ ,,F 



K+,e+ K [U K) 

\d°K\ _ _ ( 4 - 8 ) 

+ T-^(Qv(%»% ) - QK,eo(«K' U Jc)) - °- 
\ e K\ 

Indeed, we have for instance 

V K+,e+(^K( U K)) = V F K+fi+K {^ f K+e+ Y\VK{u K ))) = e +{u K ) 

We begin by introducing the following operator. For u, v € R, e° € e?°if, let 

I <9° if I 

H K ,e°(u,v) :=/4r(w) - , +| (qjf , e o («, - qjc, e o (u, n)) 
We claim that H K e o satisfies the following properties: 

—H K>e o(u,v) > 0, —H K>e o(u,v) > 0, (4.9) 

H Kje o{u,u) = n+(u). (4.10) 

The first and last properties are immediate. The second is a consequence of the CFL 
condition f)3. 12j) and the monotonicity of the method. Indeed, from the definition of 
H K e o(u, v) we may perform exactly the same calculation as in the proof of Lemma 14.21 
to prove that H K e o(u, v) is a convex combination of p^{u) and p~^(v), which in turn are 
increasing functions. This establishes the first inequality in (|4.9j) . 

We now turn to the proof of the entropy inequality (|4.8|) . Suppose first that (|4.8p is 
already established for the Kruzkov's family of entropies U(u,X) = \u — A|, F(u,X,p) = 
sgn(w — X)(f(u,p) — f(X,p)), A € K. In this case, the Kruzkov numerical entropy flux- 
functions are given by 

Qjc e o(u, v, A) := q_K,e°( u V A, i> V A) — q_K,e°( u A A, v A A). 

It is easy to check that e ° satisfies the first two conditions of the lemma. 

In fact, it is enough to prove the inequality f|4.8j) for Kruzkov's entropies only. Indeed, 
if U is a smooth function which is linear at infinity, we have (formally) 

\ [ U(u,X)U"(X)dX = \ f Tj'\u,X)U(X)dX = h5 x=u ,U(X)) = U(u), 

" il z JR z 

modulo an additive constant. Similarly, if (U, F) is a convex entropy pair, we obtain 

~ / F(u,X,p)U"(X)dX = F(u,p). 

1 JR 

Since we shall prove an L°° bound for our approximate solutions, we may suppose that 
the u above varies in a bounded set BcR. Thus, we may apply the same reasoning with 
any function which is not linear at infinity, by changing it into a linear function outside 
B. This shows that we can obtain the inequality (|4.8|) for any convex entropy pair (U, F) 
by first proving it in the special case of Kruzkov's entropies, multiplying by U"(X)/2, and 
integrating. In that case, the numerical flux will be given by 



Q K ,eo(u,v) = [ Q Ke o{u,v, X)U"{X)dX. 

1 JR 
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Again, this numerical flux satisfies the first two assumptions of the lemma, since they are 
inherited from the corresponding properties for the Kruzkov numerical flux Qj< e°( u i v > x )- 
Therefore, we now proceed to prove the inequality (|4.8p for Kruzkov 's family of en- 
tropies. This is done in two steps. First, we will show that 

_ \d°K\ 

^X+,e+( M ^' A ) ~ -|^+T"(Qx,e<'(^> n A>' A ) - Qk,A u K> u K>X)) , ] { { 

= H(u K V A, u~ K ^ V A) - H(u K A A, u K ^ Q A A). 
Second, we will see that for any u,v,X £ R, we have 

H(u V X,v V A) - H(u A A, v A A) > / K+ e+ ((/i+ v)), A) . (4.12) 

For ease of notation, we omit K, e° from the expression of H. The identity (|4.1ip and the 
inequality (|4.12|) (with u = u^, v = uZ- ) combined give (|4.8|) . for Kruzkov's entropies. 
To prove (|4.1ip . simply observe that 

= sgn(fi+(u K ) - - ^x(A)) 

= \Vk( u k) - /4( A )| 

= (^k( u k) v /4( a ) - ^(%) A VkW) 

= {^k( u k VA) -/x+(u^ A A)). 

Here, we have repeatedly used that \i\ is a monotone increasing function. The identity 
(|4.11|) now follows from the expressions of the Kruzkov numerical entropy flux, Q^ e o, 
and of H. 

Consider now the inequality (|4.12|) . We have 

H(u V A, v V A) - H(u A A, v A A) 

> (H{u, v) V H(X, A)) - (H(u, v) A H(\, A)) . 

This is a consequence of the fact that if ip is an increasing function, then (p(u V A) = 
ip(u V A) V <p(u V A) > <p(u) V <p(\), and (JOJ. Thus, we have 

H(u V A, v V A) - i?"(u AA,!)AA) 

> |ir(«,«)-£r(A,A)| 
= |fr(« lU )- A i+(A)| 

= sgn (H(u, v) - M + (A)) (tf(u, «) - /4(A)) 

= sgn (( M + J" 1 «)) - A) ( M + J" 1 (fT(u, «))) - ^ (A)) 

= ^ +e+K ((^)- 1 {H(u,v)),X). 

This establishes (|4.12p . We now choose u = WT^^v = uZ- in f|4.12[> . observe that 
H K e o [u^ , lift o ) = pL^, g0 , and combine with (|4.11|) to obtain inequality (|4.8p for Kruzkov's 
entropies. As described above, (|4.8|) will hold for all convex entropy pairs (U,F). 

To conclude the proof, observe that the functions Vk, b i n f|4.5|) are indeed convex. This 
follows from twice differentiating the representation formula 

VkM= rU'[(jM f Kie )-\u))du, (4.13) 
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and using the convexity of U and the monotonicity of (/% e ) _1 > for e = e^-,e^. This 
completes the proof of Lemma (|4.ip . □ 

4.2. Entropy dissipation estimate and L°° estimate. We now discuss the time evo- 
lution of the triangulation. As we have said, the initial hypersurface "Kq is composed of 
inflow faces e^-. We then define the hypersurfaces "K n , for n > 0, by 

and set 

X n := {K:e- K C e+ C X n }. 

An element of % n is denoted by K n when it is important to stress the time level n. 
Otherwise, when writing local in time estimates, we omit the time level n. 
Next, we introduce the following notations. For K n £ % n , write 

:= /%« = ^Cn „- > U K := u K n 

so that for instance, p^n{u^ (n ) = fi]^~ 1 (u^~ 1 ). Accordingly, we define 

VZ(fJ,):=V K (m), (4.14) 
where V^- je are the time-like entropy flux defined in Lemma 14.11 
Lemma 4.2. The finite volume approximations satisfy the L°° bound 

max \u 1 l\<( max \u° K \ + Cit n )e C2tn (4.15) 

/or some constants C\i > 0, where 

Proof. First of all, observe that from the consistency condition (|3,8p . the definition of 
/% e m (|3-2p and the divergence theorem, we have for any 

dw g f(u,p)dp = / g p (f(u,p),h(p))dp 

K n JdK n 

= |e^n|^ +1 (n) - |e^„|/i^(u) + ^ |e°|q K » e o(u, u) 

(recall that n is the interior unit normal if it is time-like, and the exterior unit normal 
if it is space-like). Moreover, with the current notations, the finite volume scheme (|3.3|) 
reads 

e o eK n 

Combining these two identities gives 

= - 7^7 / dw g f(u n K ,p)dp 

\e°\ ( 41? ) 

- Yl TO - '\(^K", e o(uK,UK e() ) -Ci K n^{u\,U n K )). 
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Next, we rewrite the right-hand side as follows: 

^ +1 « 1 ) = (i- E «^, e °)/4 +1 «)+ E ^« e o/4 +1 « e0 : 

e°£d°K n e°£d°K n 

1 div 9 f(uK,p)dp, 



(4.18) 



where 

|e°| QK"M un K> u K e0 ) - ^"^(uk^k) 



"Jf" 



This gives a convex combination of (itj^) and iJ r ^~ 1 (u 1 ^ ). Indeed, on one hand we 
have 

E a ^",e° > 

due to the monotonicity condition (|3.10p and Lemma |3,1[ On the other hand, the CFL 
condition (|3. 12j) gives 



e°eK« 



n _ n 

u K u K 



E a K n ^° < „n+l, ,.n+l, n ^ 



,n+l\-l / t :„/,,n+l\-l 



<Lip( M ^)-VLip( M ^r = 1. 
Thus, we find 



> mm K +1 (n^), min ^\u\ „)) - * / div, /(«, p)dp 
< max max J) - * / div, f(u,p)dp. 



\dj^n | J K 

Applying the monotone increasing function we find 

,n+l\-l 

n+1 ^ ■ ( n n ^ , ^fK^K 



> min min ) + Lip( ^ p — / | div 9 f{u n K ,p)\dp 
u™ +1 < max «, max < J + Lip ^ } - / | div 9 f(u n K ,p)\dp, 
rn gives 

K +1 |< max + max Lip ^f + ) /" I div„ f(u n K ,p)\dp. 



It 

v "'eOe/f" " eU/ |e 

which in turn gives 

A" 

By induction we obtain 



< max litirl+V^ max P^A' ) — /" i ^ v t (y? p)\dp. 



j=0 1 KJ 
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Now use the growth condition (|2.5p on the last term, 



V max Lip |^ ) — [ \ div g f(u J K ,p)\dp 
j^K>&0 \e+ s \ J Ki 



< v max Lip (f') V ^KCi + ^Ki) 

n 

< (Cit n + C7 2 ^ T- 7 max Iw^l). 

i=o J 

Here, the constants Ci^ may change from line to line. We have also used the fact that 

max LipQt/j!" ) _1 < C, 

which is an easy consequence of the hypotheses on the flux /. The result now follows 
from a discrete version of Gronwall inequality (see [l] Lemma 6.1]). This completes the 
proof of Lemma 14.21 □ 

Recall that if V is a convex function, then its modulus of convexity on a set S is defined 

by 

[5 := inf {V'(w) :weS}. 

Proposition 4.3. Let V£ be defined by (|4.4p . (|4.14p . and let be the modulus of 
convexity of V£ ■ Then, we have 



an+l I Oi |„+ I 

PK F \\ e K n \ l-ryn+1 n+1/ n+1 n| 2 

+ 2^ 2 \d°K n \ l^ e ° ^ ^ K j l 

; E 

E [ div g F(u" K , P ) dP+ e wr<eo 



(4.19) 



Proof. Consider the discrete entropy inequality (|4.6p . Multiplying by igo^." and sum- 
ming in if n £ X n , e° € d°K n gives 



l c C JC™ I T rn+1/— n+1 \ I + urn+l( n+1/ n \\ 

\Q0 K n\ V K W; e o)- 2^ l e ^"l^ W 

e°ed°K n 

+ E l e °l(Q^", e °«'<o) -Qx»,e° («£->«£-)) 

e g9 X n 

|e pn+1 
- |«9°if n l ^eO- 

e°ed°K n 
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Next, observe that the conservation property (|3,9p gives 



E |e°|Q/^, e o«,< ef) ) = 0. (4.21) 



K n ex n 

e°£d°K n 

Now, if V is a convex function, and if v = otjVj is a convex combination of Vj, then 
an elementary result on convex functions gives 

3 3 

Now, apply this result with the convex combination (|4.3|) and with the convex function 
V£ +1 , multiply by ffifei and sum in if n e 3C n , e° G d°.ft: n . Combine the result with 
(14301) . (OTj) to obtain 

E I4«l^ +1 « +1 « +1 ))- E 

10011.+ I 

_i_ |e HfjfnJ |-n+l _ ,,n+l^,n+l\| 2 

+ 2 '^ e0 M * 1 * j| (4.22) 

- E ie°iQ*»eo«,<)< E wr K S' 

Finally, using the identity 

div g F(u,p)dp = g p (F(u,p),h(p))dp 
JdK" 



K n JdK" 
4n\V^ n I tHu))-\e] (n \V^ n K (u))+ Y, |e°|Q^, e o(n/ 



(4.23) 



e OfZ K -n 

(with u = vOtf) yields the desired result. This completes the proof of Proposition 14.31 □ 

Corollary 4.4. Suppose that for each K € 7 h , e = e^, the function Vx, e is strictly 
convex, and that, moreover, we have 

Pk>P>0, (4-24) 

uniformly in K and n. Then one has the following global estimate for the entropy dissi- 
pation, 



N I o 1 1 + I 

EV- l e \\ e K n \ I— n+1 „n+l/„,n+l\|2 _ \ 
1^ \QO K n\ Wk^-Pk ^ U K )\ -0(*JV), (4.25) 



\d°K 

e°ed°K n 



where t^ is defined in (|4.16p . 

Proof. Summing the inequality (|4,19p for n = 0, . . . , N, we observe that the first terms 
on each side of the inequality cancel, leaving only the terms with n = and n = N. 
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Moreover, using the growth condition (|2.5p on the divergence term gives 

N on+l |„0||„+ I 

' K l e \\ e K n \ \-n+l _ ,,n+l^,n+l\| 2 

n=0K n ,e° 

< E i^oii^(^«))i+ E I4*iiv? + V£ +1 («£ +I ))i 

X°G3C /f JV + 1 eoc JV+1 



(4.26) 



JV iV i o 

+ E E i^K^i + ^KD + E E j wftf ifl S5 



n^PeX" n=0 X n G3< Tl 1 1 

The last term is estimated using (|4.2p , (|4.7p , and the growth condition (|2.5p , yielding 

E E wr<eo<E E f \x* g fWt,P)\dP 

n=0 K n eX n 1 1 n=0K n £X n K ™ 

e°£d°K n 

N 

<E E \K n \(C 1 + C 2 K\). 

n=0K n €X n 

Here, we have used that Lip V^ +1 is uniformly bounded, which is an easy consequence of 
the corresponding bounds for the flux /. The result now follows from (I4.24p and the L°° 
estimate in Lemma |4.2( which allows us to bound uniformly all the terms on the right- 
hand side of (I4.26p . Note however that this bound depends, of course, on the entropy U. 
This completes the proof of Corollary 14.41 □ 



4.3. Global entropy inequality in space and time. In this paragraph, we deduce a 
global entropy inequality from the local entropy inequality (|4,6p . This is nothing but a 
discrete version of the entropy inequality used to define a weak entropy solution. Given 
a test-function <f> defined on M rf+1 we introduce its averages 

4>{p)dp, 



le ' 



<J d°K 



e°ed°K n 



The following lemma is easily deduced from the corresponding result in the Euclidean 
space, by relying on a system of local coordinates. This result will be useful when ana- 
lyzing the approximation. 

Lemma 4.5. Let G : M — > R be a smooth function, and let e be a submanifold of M 
such that diam(e) < h. Then there exists a point p e (not necessarily in e) such that 



jG{p)dV e - G(p e ) 



<hHG\ 



We are now ready to prove the global discrete entropy inequality. 
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Proposition 4.6. Let (U, F) be a convex entropy pair, and let 4> be a non-negative test- 
function. Then, the function u h given by ()3.6[) satisfies the global entropy inequality 



oo „ 

EE/ div g {F(u%,p)<f>(p))dp 



n=0K n eX n K 

x 1 =0 



^-E E ^i^ico^r 1 ^)-^ 1 ^)) 

n=0 X n g3C" 1 1 
e°€d°K n 

00 I 0| 

+ E E T^y K ^K-€»w n+ \r£o) 

n=0 ^"gOC™ 1 

OO „ 

+ E E l(€°-<P(P))Feo(u n K , P )dp 

n=o K n eoc n e ° 

e"£d°K n 
oo „ 

"EE (/_ ttd°K-t(P))9(Hul,p),n Kne -)dp 

+ / {4> n 90 K -<j){p))g{F{u^,p),ri Kne +)dp 
+ E l ( l ) d°K9p( F (u K ,p),n K)e -)dp. 

K&OCq e K 

Proof. Prom the local entropy inequalities (|4.6p . we obtain 

\p°\\p + I / 

E l e \\^K n \ in ( (jjn+1 \ V n+l ( i, n+l (ii n \ 

K n GX n 1 1 
e°£d°K n 

+ Y |e l€o(QK»eO«,< e0 )-Qx« e o(«Sr,0) 



K n <=K n 
e°ed°K n 



l e \\ c K n \ in pn+1 

s Z> \d°K n \ K ' e °' 



Now, from the conservation property (|3.9|) we have 

E ie°ICoQ^, e o«,^ )=0 



A""g3C n 



Also, from the consistency property (|3,8p . we find 

^ |e |Q K ™ e o«,0= Yl €° [ Fe°«,p)dp 



£ / 0(p)F eO (u n K ,p)dp+ Y, [ (<K* ~ m)F e °Wk,p)dp. 



e o gc1 o X n e°ed°K n 



(4.27) 



(4.28) 
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E l e ll e K"L n v n+l(-n+l \ 
\aO K n\ (p e<> V K WK,e°) 



e°£d°K n 



Next, we have 

\d°K n \ 

E\ e \\ e K n \ in rrn+l/-n+l\ i l e \\ e k n \(±n in u m+l/-n+l\ 

\Q0 K n\ ^kVk fc£ e o) + ^ [flO^nl ~ ^kWk (/%>) 

Here, we have used the fact that for a convex function V and a convex combination 
v = Y^j o/.jVj one has 

v(v) <e«^)> 

with the convex function V^3 +1 and the convex combination (14. 3p . Also, 

Therefore, the inequality (|4.28|) becomes 
/Per 

K n £X n e ° 
e°ed°K n 

* E wrCo^- e wr^o-fc)^ 1 ^) (4-29) 



+ E [ o (€o-Hp))F e o(u n K ,p)dp 



gOggO^n 



=: A h + B h + C h 
The first term in (|4.29p can be written as 

E / + ^0») 9{F{u n K +1 , p) - F« , p) , n K „ e + (p)) dp 
+ E / + (^-^(p))5(FK + \p)-F(<,p),ri i ,„ i g+(p))dp. 
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Combining this result with the identity 



div g (F(u,p)(p{p)) dp = / 4>(p) g(F(u,p),n aK )dp 

JdK 



(f>(p)g{F{u,p),n Ke +)dp+ <p{p) g(F(u,p),n Ke -)dp 

Je K 

+ ^ 4>(p) F e°( u ,P)dp 



e°ed°K 

(with u = u K ) and in view of H4.29H we see that 
div g (F{u n K ,p)cf>{p))dp 
<A h + B h + C h 



( / M9{n^ 1 ,p),n K n ie+K (p))dp+ ! cf>(p)g(F(u, P ),n Ke -)d P ^ 

E / + ^K-m)g{nu K + \p)-F{u n Kl p) 1 i, K ^ e+K {p))dp. 



The inequality (|4.27p is now obtained by summation in n. First, the (summed) terms 
A,B,C give the first three terms on the right-hand side of (|4.27p . and by discrete inte- 
gration by parts in the last terms of the above inequality, we find the remaining terms. 
This completes the proof of Proposition 14.61 □ 



5. Proof of convergence 

This section contains a proof of the convergence of the finite volume method, and is 
based on the framework of measure-valued solutions to conservation laws, introduced 
by DiPerna [1] and extended to manifolds by Ben-Artzi and LeFloch [2]. The basic 
strategy will be to rely on the discrete entropy inequality (|4.27[) as well as on the entropy 
dissipation estimate (|4.25p . in order to check that any Young measure associated with 
the approximate solution is a measure-valued solution to the Cauchy problem under 
consideration. In turn, by the uniqueness result for measure-valued solutions it follows 
that, in fact, this solution is the unique weak entropy solution of the problem under 
consideration. 

Since the sequence u h is uniformly bounded in L°°(M), we can associate a subsequence 
and a Young measure v : M — > Prob(M), which is a family of probability measures in R 
parametrized by p G M. The Young measure allows us to determine all weak-* limits 
of composite functions a(u h ), for arbitrary real continuous functions o, according to the 
following property : 

a(u h ) (u,a) as h — > (5.1) 

where we use the notation 

(u, a) := / a(X)dv(X). 
Jr 

In view of the above property, the passage to the limit in the left-hand side of (|4.27p is 
(almost) immediate. The uniqueness theorem [H [2] tells us that once we know that v is 
a measure- valued solution to the conservation law, we can prove that the support of each 
probability measure v p actually reduces to a single value u(p), if the same is true on IHo, 
that is, Vp is the Dirac measure 5 u ( p y It is then standard to deduce that the convergence 
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in (|5.ip is actually strong, and that, in particular, u h converges strongly to u which in 
turn is the unique entropy solution of the Cauchy problem under consideration. 

Lemma 5.1. Let v v be the Young measure associated with the sequence u h . Then, for 
every convex entropy pair (U, F) and every non-negative test-function <f> defined on M 
with compact support, we have 



(v p ,dbr g F{;p))<j)(p) +g{(v P ,F(; P )), V0) dp 



4>(p)g{{v p ,F(-,p)),nx )dp+ 4>(p){v P ,U'(-)div g f(-,p)}dp<0. (5.2) 



M 



Proof of Theorem \3.4\ According to the inequality (15. 2D , we have for all convex entropy 
pairs (U, F), 

div g (v,F(-)) - (u, (div 9 F)(-)) + (v,U'(-)(div g /)(■)) < 

in the sense of distributions in M. Since on the initial hypersurface "Kq, the (trace of 
the) Young measure v is the Dirac mass 5 U0 (because no is a bounded function), from 
the theory in [2] there exists a unique function u E L°°(M) such that the measure v 
remains the Dirac mass 5 U for all Cauchy hypersurfaces "K^ < t < T. Moreover, this 
implies that the approximations u h converge strongly to u on compact sets at least. This 
concludes the proof. □ 

Proof of Lemma \5.1l The proof consists of passing to the limit the inequality (|4.27p and 
using the property (|5.ip of the Young measure. First, note that the left-hand side of the 
inequality (j4.27j) converges immediately to the first integral term of (|5.2p , in view of (|5.ip . 
Next, take the last term of (I4.27|) . Using again (|5.ip and the fact that <Pqo k — <t>{p) = 
(t + h) , we see that this term converges to the second integral term in (|5.2p . 

Next, we will prove that the first term on the right-hand side of (|4.27p converges to 
the last term in (15.21). Observe first that 



\ e K n I J K n 

Therefore, we obtain 

00 I 0| 

-EE ^i^ico^r 1 ^)-^ 1 ^)) 

n=0 K n £X n 1 1 
e°ed°K n 

oo I o I r- 

= E E WH^o(W +1 (Ateo)/ div g f{ul,p)dp) 

n =o x»ex" 1 1 Rn 

e°ed°K n 



+ \e + Kn \0{r 2 \ 



oo i o I r- 

E E |^^((W +1 (Ateo)-^K +1 )) / div, 



\e + Kn \0(r 2 ) + U'{u n K +1 ) / dW g f(u n K ,p)dp)) 



Now, note that from the expression of V, (|4.13p 



< supU" sup Lip^)" 1 !^ - M «+ 

n,K n 
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and so, using the L°° bound (|4.15p and the growth condition (|2.5|) . we find 

co I 0| r 

E E |^Co(W +1 (^ e o)-^(< +1 )) ltiv g f(u n K,p)dp 

e ° e9 °^' 1 (5.3) 

< —±lJ—rh n \ K n \\u n+1 — li n+1 \ 

n=0 X^GOC 1 1 1 

Applying Cauchy-Schwartz inequality and the entropy dissipation estimate (|4.25p . we 
find that this term tends to zero with h. We are left with the term 

CO . 

E E ^kU'^k 1 ) / div g f(u n K ,p)dp), 
which is easily seen to be of the form 

CO ~ ~ 

E E U'( U n K ) J( P )dw g f(u n K ,p)dp + 0(h)^ (v p ,U'(-)div g f(;p))d P . 

It remains to check that the remaining terms in (|4.27p tend to zero with h. Namely, 
the second term on the right-hand side can be written as 

00 I 0j 

E E T^\ e U(<(> n d OK-€o)V n+ \TI^o) 

n=0 K n £X n 1 1 
e°ed K n 

00 (5 4) 

= E E i^i^Kfc-^)^ 1 ^)-^^^ 1 ^ 1 ))) : ; 

n=0 K n £X" 1 1 
e°€d°K n 

= o(l), 

by Cauchy-Schwartz inequality and the entropy dissipation estimate (|4,25p . Next, the 
second term on the right-hand side of (|4.27p satisfies 



OO „ 

E E m-Hp))F e o(u K , P )dp 

n=0 K n £3C n 
e°ed°K" 

oo „ „ 

E E (€o-<P(p))(Fe°(u n K,p)-} F e o(u n K ,q)dq)dp. 

n Tsn^cvn Je° Je° 



ra =o rer 

e°£d°K n 

In view of the regularity of 4> and F, this term is bounded by 

CO 

E \d°K\(D(r K n + h) 2 . 

Using the CFL condition (|3.12p . and the property (|3.1ip . we can further bound this term 
by 

CO 

E E \ein\°( T K"){0(T K "+h) + 0(h 2 /T K n)} = o(l). 

n=0K n €X n 
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Only the term 

A*):=E E ( [_ (<P n d ^-<P(p))9(F(u n K ,p),n Kn eK )dp 

n=0 K n £X n \ Je K™ 

+ f + (<t>&K-<t>(p))9{F{ul,p)^ Kne+K )dp\ (5.6) 

remains to be bounded. Here, we will use that our triangulation is admissible, in the 
sense of Definition 13.31 First of all, observe that using Lemma [4. 51 one may replace <f>go K 
by 4>{p°Kri)i and 4>qok ^ 4>(p°k«-i)i where p° Kj denotes the center of c^iP, with an error 
term of the form C(r + /i)||^>||c2 ||i^||i,<». Next, we replace 

with 

HPK)9(F{u^,p+),n Kne +{p+)), 

and similarly for e^-. Using the property (I3.15p . the resulting error term is seen to be of 
the form C(r + h) \\4>\\c 2 ll-^llc 2 ■ We then have 



+ |e^|(<KPA'"-0 - < / , (P^n-i))5'(^('"^,^ n - 1 ),n^ e - ( (p+ n _ 1 ))^ 

+ (T + h)\\<j>\\ (P {\\F\\ L - + \\F\\ C 2). 

Now, performing a Taylor expansion of <j> and using the definition of wk (recall that wk 
is the future-oriented unit vector at p~^- tangent to the geodesic connecting and p® K ) 
we find, for instance, 

Hpk) ~ <P(Pk) < hg(w K ,V<p(p+)) + 0(h 2 ). 

Therefore, by the definition of E(K) and using (|3.13p . we may express this using the 
Cartesian deviation of the triangulation, 

CO / 

A0)<E E h l4l5(w^,V^(p+)) 5 (F«,p+),n^ 4 (p+„)) 

n=0if»g3C" V 

+ l e A'kK»-!> V0(pJ„_ 1 ))3(-F(uSr,pJ n _ 1 ),n ifn e - (pj„_0)l 

+ (T + h)\\<t>\\ c2 (\\F\\ L ~ + \\F\\ c2 ) 
<fc E {\K\e,(K)-\K-\Z(K-))(V<t>,F(u h )) 

K n £l h 

+ {T+h)U\\ C 2{\\F\\ L ~ + \\F\\ C 2) 

< ri(h)\\<j>\\ci\\F\\v° + (r + h)U\\ C 2(\\F\\ L ~ + ||F|| C 2), 
which tends to zero since Tj(h) — > 0. This completes the proof of Lemma l5.11 □ 
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6. Examples and remarks 

6.1. A geometric condition. The following proposition gives a condition of a geometric 
nature which is sufficient for a triangulation to be admissible in the sense of ()3.13j) . Recall 
that denotes the center of mass of e^-, and that the vector w denotes the tangent at 
e\, to the geodesic from e^- to the center of d°K. 

Proposition 6.1. Let7 h be a triangulation. Suppose that, for each element K, the scaled 
exterior normals \e^\r\. K e + and \e~^\n K e - , and the vectors Wx andw^- approach in the 
limit, in the following weak sense: for every smooth vector field X , one has 

n(h) 

\9(,\e K \n K>e pX) - g(\e K \n Ke -,X)\ < ^ |tf |pr|| £oo(J0 , (6.1) 

\g(w K ,X) - g(w K -,X)\ <^t k \\X\\ Loo(k) , (6.2) 

where the expressions are evaluated at the centers of mass of e^ and e~^-, and r](h) is such 
that n{h) — > 0. Then, T h is an admissible triangulation in the sense of (13. 13|) . 

Proof. Let <E> be a smooth vector field, and let ^ be smooth on each element K. We have 

{\K\E(K)-\K-\E{K-))(*,* /K ) 

K£l h 

= Yl (I4l w ^® n i^ e + - \ e ~ K \™K- ® n K, £ - 
K£7 h 

= Yl f( w ^' $ )f(|e^|n Ke +,* /if ) -g(w K -,$)g(\e- K \n Ke -,-$ /K ) 
= 9(^K-,^)(9(\^\n Ke p^ /K )-g(\e],\n Ke -,^ /K )) 

K£l h 

+ (s(wjf,$) - g(w K -,^)) g{\e^\n Ke +,^ /K ) 

<r]{h)/h IKWMl^W^ /k\\l~> <ri(h)/h\\$\\ L ™ sup \\^> /K \\ L ™. 
KeT h K 

In view of (|3.13j) . this shows that 7 h is an admissible triangulation. □ 

For instance, one can easily check that if a triangulation is subordinate to a given 
foliation (in the sense that each connected component of the set of all outgoing faces 

{e+:Ke 7 h } 

is contained in a certain Cauchy surface), and if, moreover, each lateral face e° is ev- 
erywhere tangent to a given, fixed, smooth time-like vector field, then the hypotheses of 
Proposition 16. II hold. However, our conditions (|3.13p or (|6.1[) . (|6.2p allow for more general 
triangulations, which need not satisfy such regularity assumptions. 

6.2. Lax Friedrichs-type flux-functions. Our general space-time setting requires some 
care when defining particular numerical fluxes. In particular, the absence of a canonically 
singled out conservative variable implies that one needs to use the values of ^ Ke {u), for 
both K and its neighboring element K e o, to compute the flux along an interface. For 
instance, consider the following natural generalization of the Lax-Friedrichs-type flux, 

ciK, e o(u,v) = -(^ Ke0 (u)+^ Ke0 (v)) + —^(fM f + (u)-^ + {v)), (6.3) 

Z ' Z 5 K e'J ' K n 
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where the constants D K e o satisfy 

le+l 

This numerical flux is conservative, and it is monotone, as may be checked using the CFL 
condition f)3. 12j) . The consistency property (|3.8p is seen to be valid once we remark that 
this numerical flux can be written as 

c\ Kte p(u,v,m(u),H2(v)) = ^(/it ie p(«)+A*t ie p(u)) + — Tp-(A*i(«) -yU 2 (V)), 
and that the consistency property reads, in fact, 

<\ K ^(u,u,ii{u),ii{u)) = H f K ^{u). 

Acknowledgements 

The authors were partially supported by the Agence Nationale de la Recherche (A.N.R.) 
through the grant 06-2-134423 entitled "Mathematical Methods in General Relativity" 
(MATH-GR), and by the Centre National de la Recherche Scientifique (CNRS). The first 
author (P.A.) was also supported by the FCT-Fundaqao para a Ciencia e Tecnologia (Por- 
tuguese Foundation for Science and Technology) through the grant SFRH/BD/17271/2004. 

References 

[1] Amorim P., Ben-Artzi M. and LeFloch P.G., Hyperbolic conservation laws on manifolds: total 
variation estimates and the Bnite volume method, Meth. Appl. Anal. 12 (2005), 291-324 

[2] Ben-Artzi M. and LeFloch P.G., The well-posedness theory for geometry compatible, hyperbolic 
conservation laws on manifolds, Ann. Inst. H. Poincare: Nonlin. Anal. (2007). 

[3] Cockburn B., Coquel F., and LeFloch P.G., Convergence of Enite volume methods for multi- 
dimensional conservation laws, SIAM J. Numer. Anal. 32 (1995), 687-705. 

[4] DiPerna R. J., Measure- valued solutions to conservation laws, Arch. Rational Mech. Anal. 88 (1985), 
223-270. 

[5] Kruzkov S.N., First-order quasilinear equations with several space variables, Math. USSR Sb. 10 
(1970), 217-243. 

[6] LeFloch P.G., Neves W., and Okutmustur B., in preparation. 
[7] LeFloch P.G. and Okutmustur B., in preparation. 

[8] Panov E.Y., On the Cauchy problem for a first-order quasilinear equation on a manifold, Differential 
Equations 33 (1997), 257-266. 

P. Amorim, P.G. LeFloch, and B. Okutmustur, Laboratoire Jacques-Louis Lions & Centre 
National de la Recherche Scientifique, Universite de Paris 6, 4 Place Jussieu, 75252 Paris, 
France. 

E-MAIL : Amorim@ann.jussieu.fr, LeFloch@ann.jussieu.fr, Okutmustur@ann.jussieu.fr. 



